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We show that the so-called flat-space rotational Killing vector method for measuring the Cartesian 
components of a black hole spin can be derived from the surface integral of Weinberg’s pseudoten¬ 
sor over the apparent horizon surface when using Gaussian normal coordinates in the integration. 
Moreover, the integration of the pseudotensor in this gauge yields the Komar angular momentum 
integral in a foliation adapted to the axisymmetry of the spacetime. As a result, the method does 
not explicitly depend on the evolved lapse a and shift /3* on the respective timeslice, as they are 
fixed to Gaussian normal coordinates, while leaving the coordinate labels of the spatial metric 7 ij 
and the extrinsic curvature Kij unchanged. Such gauge fixing endows the method with coordi¬ 
nate invariance, which is not present in integral expressions using Weinberg’s pseudotensor, as they 
normally rely on the explicit use of Cartesian coordinates. 
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I. INTRODUCTION 


Since the binary black hole (BBH) merger break¬ 
through simulations of about a decade ago [lH3|, ever¬ 
growing computational resources and advances in the 
numerical methods used to simulate these systems have 
made the exploration of the vast initial parameter space 
possible (see e.g. Q and references therein for a recent 
overview of the status of BBH simulations). The initial 
parameters of BBH simulations are the BH mass ratio 
and the six components of their initial spin vectors. The 
investigation of these initial parameters has led to signif¬ 
icant discoveries, as the occurrence of the orbital hang¬ 
up i and the presence of the so-called super-kicks where 
the final BH is displaced from the orbital plane after its 
formation Q. 

The BH spin and in particular its orientation may also 
play a non-negligible role in non-vacuum spacetimes in¬ 
volving BHs surrounded by matter as, e.g. in the form of 
accretion disks, a situation commonly encountered in bi¬ 
nary neutron star merger simulations. Recently, we have 
performed numerical relativity simulations of tilted self- 
gravitating accretion disks around BHs, investigating the 
precession and nutation the BH undergoes as it accretes 
mass and angular momentum from the torus Q. In order 
to carry out a quantitative analysis of such kind of sim¬ 
ulations, it is obviously necessary to measure both the 
magnitude of the BH spin and also its direction in space. 
Such a study inspired the work we present in this paper. 

One of the standard methods in numerical relativity to 
measure the magnitude of the angular momentum of the 
BH horizon is described in Q. This method is based on 
the so-called isolated horizon formalism Q and the gen¬ 
eralization to dynamical horizons M- In this approach 
the BH spin is calculated by performing the following 








Jah = ^ ^ {^R^Kab) dS, (1) 

where '0“ is an approximate rotational Killing vector on 
the horizon that has to be determined numerically (see 
0 for a method of finding ■)/)“), is the outward point¬ 
ing unit vector normal to the horizon, Kab is the extrinsic 
curvature on the horizon surface and dS is the surface ele¬ 
ment. This method does not, however, give the direction 
of the BH spin in the 3D Cartesian reference frame of the 
computational grid of a numerical relativity simulation. 

The direction of the BH spin in numerical relativity is 
commonly measured by the approach suggested by [111. 
In this approach the BH spin direction is simply de¬ 
fined as the Euclidean unit vector tangent to the coor¬ 
dinate line joining the two poles on the horizon (i.e. the 
two points where the axially symmetric vector ■0“ van¬ 
ishes). The approximate Killing vector field 0“ on the 
horizon is obtained numerically using spherical-polar co¬ 
ordinates and the accuracy in the spin direction is typi¬ 
cally about a few angular grid zones. This definition of 
the spin vector reproduces the Bowen-York spin param¬ 
eter on the initial slice and gives satisfactory results as 
long as the BH horizon does not become too distorted. 
Moreover, [HI present another method for finding the 
spin magnitude and direction, using flat-space coordi¬ 
nate rotational Killing vectors to calculate the Cartesian 
components of the BH spin and its magnitude from the 
Euclidean norm of the resulting vector. The flat-space 
Killing vector method has the practical advantage that 
the vector 0“ used in the surface integral o is given an¬ 
alytically and is constant, therefore it does not have to 
be found numerically on each timeslice. 

In this paper we show how the flat-space Killing vector 
method can be derived by performing a surface integral of 
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Weinberg’s energy-momentum pseudotensor [l^ . By us¬ 
ing the 3-1-1 split of spacetime and Gaussian coordinates, 
it is possible to express the angular momentum of a given 
volume using Weinberg’s energy-momentum pseudoten¬ 
sor in a simple form that allows for a straightforward 
calculation of the spin vector of the BH horizon. Wein¬ 
berg’s energy-momentum pseudotensor is a symmetric 
pseudotensor derived by writing Einstein’s equations 
using a coordinate system that is quasi-Minkowskian, 
i.e. with the four-dimensional metric g^i, approaching the 
Minkowski metric at infinity. Although it is not gen¬ 
erally covariant, the pseudotensor is Lorentz covariant, 
and with the appropriate choice of coordinates it pro¬ 
vides a measure of the total angular momentum of the 
system. In the following Greek indices run from 0 to 3 
while Latin indices run from 1 to 3. We use geometrized 
units (G = c = 1) throughout. 

II. 3+1 SURFACE INTEGRAL OF 
WEINBERGS’S PSEUDOTENSOR IN GAUSSIAN 
COORDINATES 

In this section, we first briefly review the calculation 
of the angular momentum contained in a volume using 
Weinberg’s pseudotensor. Next, we express the resulting 
surface integral in terms of the 3+1 spacetime variables 
on a given timeslice. Finally, we show that by choosing 
Gaussian coordinates the integral reduces in complexity 
and is analytically equivalent to the flat-space rotational 
Killing vector method. 


and indices of linearized quantities are raised and lowered 
with ? 7 ^y. 

Using the pseudotensor, the volume integrals giving 
the total four-momentum of the volume are given by 



Furthermore, the pseudotensor defined by Eq. @ 
is symmetric, which allows one to use it to calculate the 
total angular momentum in a volume V using the follow¬ 
ing volume integral: 


= 


Jv 


IV 

1 

Stt 


- X a x. 

dx'- dx'- J 


( 6 ) 


As Weinberg remarks, the physically interesting Garte- 
sian components of the angular momentum contained in 
the volume are 
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Using Gauss’ law the volume integral can be trans¬ 
formed to the following surface integral over the bound¬ 
ing surface: 
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A. Angular momentum with Weinberg’s 
pseudotensor 


Weinberg’s energy-momentum pseudotensor is ob¬ 
tained by writing the Einstein equations in a coordinate 
system that is quasi-Minkowskian in Cartesian coordi¬ 
nates, so that the metric approaches the Cartesian 
Minkowski metric ? 7 ^j/ = diag(—1,1,1,1) at infinity as 
follows 


9tiu = 




( 2 ) 


where hfiy does not necessarily have to be small every¬ 
where. Then, by writing the Einstein equations in parts 
linear in one arrives at an energy-momentum pseu¬ 
dotensor , which is the total energy-momentum “ten¬ 
sor” of the matter fields, Tx,^, and of the gravitational 
field, Gk, 

[Tx. + tx.] = ( 3 ) 

where is the superpotential given by 
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where rd is the unit normal to the surface of integration 
and dS the surface element. 

The convergence of the four-momentum volume inte¬ 
grals involving the pseudotensor critically de¬ 
pends on the rate at which the metric g^^, approaches 
the Minkowski reference metric at large distances. Given 
the following behaviour of as r ^ oo. 
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where r = (x^ + y'^ + it can be shown that the 

energy-momentum “tensor” of the gravitational field, 
behaves at large distances as 


V = 0{r-^), (10) 


which in turn shows that the four-momentum volume 
integral (0 converges. The convergence of the total an¬ 
gular momentum volume integral ([6]) and of the corre¬ 
sponding surface integral ([51) is more problematic, due to 
the appearance of x^ in the volume integral. This is also 
observed in the convergence properties of the integrals 
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of the ADM quantities [l^, where the surface integrals 
for the ADM mass and linear momentum converge when 
imposing fall-off conditions like those of Eq. while 
the calculation of the ADM angular momentum gener¬ 
ally requires stronger asymptotic fall-off conditions M- 
We shall return to the issue of the convergence of Eq. ([51) 
after we have expressed it in terms of the 3-1-1 variables 
and in Gaussian normal coordinates in the next section. 


B. The angular momentum pseudotensor integral 
in Gaussian coordinates 

We can express the total angular momentum given by 
Eq. ([ 8 ]) in Gaussian normal coordinates (also called syn¬ 
chronous coordinates), which represent free-falling ob¬ 
servers. We start by doing a 3-1-1 decomposition of the 
four-dimensional metric 
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lijP^ 
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( 11 ) 


where a is the lapse function, /3® the shift vector, and 7 ^ 
the spatial metric induced on the hypersurface. From the 
requirement that the metric approaches Gartesian 
Minkowski space at infinity ®, we see that 
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If we now express the angular momentum surface inte¬ 
gral, Eq. (| 8 |), in terms of the 3-1-1 variables we find that 
can be written as 


= - 


1 

IGtt 


djljmPn , 

dx^ dx^ 


^ d{ljk - Sjk) ^ di'-fik - Sik) 


dt 


dt 


+ ijmrSu - n^dS. (13) 


Moreover, in terms of the 3-1-1 variables, Gaussian co¬ 
ordinates are defined by the following choice of the lapse 
and shift vector: 


a = l, (3^=0, (14) 

so that Hqq = hoi = hio = 0. In this gauge, Eq. (fl^ 
considerably simplifies to 


where Cp is the Lie derivative with respect to the shift 
vector f3\ to see that the time derivative of the spatial 
metric d^ij/dt in Gaussian coordinates is simply 


dlij 

dt 


- 2 iLy. 


(17) 


Substituting Eq. (fTTll in Eq. (fT51) . we find that 


^ JJ - XjKik)n'^dS. (18) 

Finally, using Eq. 0 , the three components of the 
Gartesian angular momentum vector of a volume are 
given by 

J. = ^ iv^sk - zK^k) n'^dS , 

Jy = J^^ = -^JJ {zKik - xK3k)n^dS, (19) 

= = i^K2k - yKik) n^dS . 

Introducing the components of the three Gartesian 
Killing vectors of the rotational symmetry of Minkowski 
space 


= ( 0 ,-z,y) 

Cy = -x) ( 20 ) 

= [-y,x,Q) 

we can rewrite the surface integrals of the three Gartesian 
components of the angular momentum in the following 
way: 


J^ = ^JJ^K,k{^,yn'^dS. ( 21 ) 

Thus, Weinberg’s identification of the (2,3), (3,1) and 
( 1 , 2 ) components as being the physically interesting ones 
is now clearly seen from Eq. (ED), as it is the rotational 
Killing vectors of Minkowski space that enter in the cal¬ 
culation of the Gartesian components of the total angular 
momentum of the volume. 

Note that this form of the angular momentum is re¬ 
markably similar to that of the ADM angular momentum 

0 : 


Ji = — lim 

StT r-^-oo 



{K,k - K^jk) {^iYu^dS . 


( 22 ) 




We can now use the definition of the extrinsic curvature 


K, 
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If the integration is done over a sphere, the components 
of the surface normal are given by 

(23) 

\r r r/ 

so that (^i)-^ and are orthogonal vectors, 

loki^iYu^ = {£,i)^nk =0 Vi. (24) 
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Therefore, the part of the integral containing the trace of 
Kij in Eq. (IMl) vanishes for spherical surfaces and there¬ 
fore equations (EH) and (1^ are identical. We have thus 
shown that by using Weinberg’s pseudotensor in Gaus¬ 
sian coordinates we obtain the total ADM angular mo¬ 
mentum evaluated at spatial infinity, when the integra¬ 
tion surface is a sphere. We might still need to impose 
a stricter asymptotic behaviour than the asymptotic Eu¬ 
clidean flatness of [l^ (for instance the quasi-isotropic 
or asymptotic maximal gauge), but as [I^ noted, the 
Kjki^iYn^ part of Eq. (1^ converges in practice. We 
are, however, interested in evaluating Eq. m quasi- 
locally, that is, associated with finite 2-surfaces (in our 
actual applications, these will be apparent horizons of 
black holes 0 )- 

For an axisymmetric spacetime, the angular momen¬ 
tum can be calculated via the so-called Komar angular 
momentum Il6ll, w hich is defined as (following again the 
notation of |14.flH|l: 

= (25) 

where cj)" is the axial Killing vector. Note the extra factor 
of 2 in the denominator, known as Komar’s anomalous 
factor 0. The Komar angular momentum integral does 
not have to be evaluated at spatial infinity, but is valid 
for every surface. In [l3, [13 it is shown that using a 
slicing adapted to the axisymmetry of the spacetime, and 
expressing Eq. EH in terms of the 3-1-1 variables, the 
Komar angular momentum becomes 

Jk = ^ (26) 

In [13 the above integral is evaluated for a Kerr BH 
in spherical Boyer-Lindquist coordinates, and the angu¬ 
lar momentum is found to be Jk = Ma, as expected, 
where M and a are the black hole mass and spin pa¬ 
rameter, respectively. As the two integrals (EH and (IMl) 
have exactly the same structure, and the latter is co¬ 
ordinate (but not foliation) invariant, we arrive at the 
conclusion that the introduction of Gaussian coordinates 
has led to a coordinate invariant expression for the an¬ 
gular momentum derived from Weinberg’s pseudotensor, 
namely the Komar angular momentum. Note the ab¬ 
sence of the anomalous factor of 2 in our final expression 
(EH- It therefore seems that it is possible to relax the 
restriction of using Gartesian coordinates in calculations 
involving Weinberg’s pseudotensor. 

C. Measuring the angular momentum in numerical 
relativity simulations 

It is easy to check that not only the choice of Gaussian 
coordinates simplifies the calculation of the total angular 
momentum via Weinberg’s pseudotensor, but also that it 
makes straightforward the implementation of the above 


expressions in a numerical relativity 3D Gartesian code 
based on the 3-1-1 decomposition. For instance, if using 
the widely adopted BSSN formulation [isi - l^ , the extrin¬ 
sic curvature Kij of the spatial slices is closely related to 
one of the evolved variables, namely the traceless part 
of the conformally related extrinsic curvature. We note 
that in present-day numerical relativity simulations one 
does not typically use Gaussian coordinates for the ac¬ 
tual numerical evolutions. This has to do with the fact 
that Gaussian coordinates can only be used in the close 
vicinity of a spatial hypersurface, as the geodesics ema¬ 
nating from the hypersurfaces will eventually cross and 
form caustics in a finite time [13 • Furthermore the foli¬ 
ation is not singularity-avoiding, which means Gaussian 
coordinates are unsuitable for the numerical evolution of 
spacetimes containing curvature singularities. Instead, 
the gauge conditions most commonly employed today in 
numerical relativity belong to the family of the so-called 
moving puncture gauges, which consist of the “l-|-log” 
condition for the lapse function E^and the Gamma 
driver condition for the shift vector [22|| . However, one 
can use the numerical solution for the extrinsic curvature 
Kij in Eq. (EH due to the freedom to choose any gauge 
for calculations done on each timeslice. 

In addition, Eq. (1^ is actually equivalent to the 
method proposed by [li| for the calculation of the an¬ 
gular momentum of a volume using flat-space coordinate 
rotational Killing vectors (cf. Eq. (EH)- To see this, con¬ 
sider the definition of the Killing vectors in Cartesian 
coordinates given by [TTj |: 

'^x = Zc),{y-yc)], 

V’y = [(^-^c),0,-(x-Xc)], (27) 

i’t = [-{y -yc),{x -Xc),% 

where (Xc, yc, Zc) is the coordinate centroid of the appar¬ 
ent horizon, which has to be subtracted to avoid including 
contributions from a possible orbital angular momentum 
of the BH about the center of the computational grid in 
the calculation of its spin. Upon substituting their flat- 
space coordinate rotational Killing vectors into Eq. o, 
we find that 

^ JJ - zK2b) n^dS , 

^ JJ - xKsb) n^dS , (28) 

^ JJ i.xK2b - yKib) n^’dS , 

where we have set Xc = yc = Zc = 0 for simplicity. We see 
that the two sets of expressions for the Cartesian compo¬ 
nents of the angular momentum vector of the AH, those 
from Weinberg’s pseudotensor evaluated in Gaussian co¬ 
ordinates and those from the flat space rotational Killing 
vector method, are equivalent and equal to the Komar 
angular momentum in an axisymmetric spacetime. 
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III. DISCUSSION 


As we have shown, the flat-space rotational Killing vec¬ 
tor method of [ll| can be derived from Weinberg’s pseu¬ 
dotensor when using Gaussian coordinates. These coor¬ 
dinates have two interesting properties that make them 
particularly useful for the evaluation of the angular mo¬ 
mentum pseudotensor integral, Eq. (I8|). First, as we have 
shown, the complicated integral ([5]) reduces to the much 
simpler expressions given by Eq. m and this final ex¬ 
pression is equal to the Komar angular momentum in¬ 
tegral in a foliation adapted to the axisymmetry of the 
system. As a result, one does not need the knowledge of 
the shift vector and of its spatial derivatives on the sur¬ 
face of integration, which in practice would involve more 
quantities that one would need to interpolate onto the 
horizon surface for the calculation of the spin, thus also 
avoiding the numerical error associated with the com¬ 
putation of the finite difference approximation to those 
spatial derivatives. Second, Gaussian coordinates triv¬ 
ially satisfy the necessary falloff conditions for the lapse 
and shift. Moreover, by using Gaussian coordinates we 
recover the ADM angular momentum evaluated at spa¬ 
tial infinity, provided we use a spherical surface of inte¬ 
gration. 

It is generally known that the various energy- 
momentum pseudotensors proposed in the literature are 
not covariant and care has to be taken when evaluat¬ 
ing them in different coordinate systems and gauges. 
(See [ 2 ^ for a review on quasi-local mass and angular 
momentum in General Relativity, where the problems 
arising when using pseudotensors for the calculation of 
mass and angular momentum are also discussed.) The 
derivation of Weinberg’s pseudotensor relies crucially on 
the reference space being Cartesian Minkowski. In his 
textbook [I^ Weinberg states that a spherical polar co¬ 
ordinate system would lead to a gravitational energy den¬ 
sity concentrated at infinity. While being non-covariant 
is generally not desirable, Weinberg’s method is em¬ 
ployed in a Cartesian grid, and Gaussian coordinates 
guarantee the correct asymptotic behaviour of the lapse 
and shift, irrespective of the asymptotic behaviour the 
evolved lapse and shift may posses, which are, as pre¬ 
viously stated, not explicitly used in the calculation of 
the AH spin on the respective timeslice. Furthermore, 
we have shown that Gaussian coordinates transform the 
pseudotensor angular momentum surface integral (I13II to 
the Komar angular momentum integral (1261) which is co¬ 
ordinate independent. The use of Gaussian coordinates 
(as an explicit gauge-fixing) seems to therefore remove 
the coordinate restrictions of the pseudotensor. 

When ll| introduced the flat-space rotational Killing 
vectors for the calculation of the BH spin direction, the 
authors stated that they could not guarantee the correct 


results for all times because the method is not gauge in¬ 
variant. However, as we have seen, such method can be 
derived from the integration of Weinberg’s total angu¬ 
lar momentum pseudotensor over the apparent horizon 
surface when using Gaussian normal coordinates in the 
integration. As a result, the method does not depend on 
the evolving lapse and shift, as the gauge is fixed to the 
Gaussian normal coordinates on the respective timeslice. 
We stress that the evolution of the lapse and shift during 
the free evolution of the spacetime does not enter the cal¬ 
culation, given the coordinates evolve in such a way that 
an AH is found at all times during the evolution, which 
is usually the case in puncture evolutions with the BSSN 
system. There is a dependence on the gauge evolution 
via the extrinsic curvature Kij that is interpolated onto 
the AH for the calculation of the spin direction, but the 
same is true for the expression of the spin magnitude in 
Eq. 

In the authors have shown that Eq. ([1]) will give the 
spin magnitude provided an approximate Killing vector 
can be found on the horizon and is gauge independent on 
the respective time-slice if the approximate Killing vec¬ 
tor field '0“ is divergence-free. Here we have shown that 
both methods (i.e. either via Weinberg’s pseudotensor in 
Gaussian coordinates or via flat-space rotational Killing 
vectors) yield the Komar angular momentum when the 
latter is expressed in a foliation adapted to the axisym¬ 
metry. We note that the restriction to axisymmetry turns 
out in practice not to be a major weakness, as numerical 
relativity simulations repeatedly show that the remnants 
of binary black hole mergers and perturbed Kerr black 
holes typically settle down to the axisymmetric Kerr so¬ 
lution quickly [ 23 , 123 . Moreover, both methods provide 
a measure of the BH spin magnitude and direction that 
is not explicitly dependent on the lapse and the shift on 
the respective time-slice. As both methods use the fixed 
rotational Killing vectors of Minkowski space, they mea¬ 
sure the spin contribution from the axisymmetry of the 
AH. 
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